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ABSTRACT 

Blandford-Znajek (BZ) mechanism has usually been studied in the literature for accre- 
tion with considerably high angular momentum leading to the formation of a Keplerian 
disc. However, in nearby elliptical galaxies, as well as for our own Galactic centre, low 
angular momentum accretion is prevalent. Such sub-Keplerian accretion has complex 
dynamical features and can accommodate stationary shocks. In this letter, we present 
our calculation for the maximum efficiency obtainable through the BZ mechanism for 
complete general relativistic low angular momentum flow in the Kerr metric. Both 
shocked and shock free flow has been studied in detail for rotating and counter rotat- 
ing accretion. Such study has never been done in the literature before. We find that 
the energy extraction efficiency is low, about 0.1%, and increases by a factor 15 if the 
ram pressure is included. Such an efficiency is still much higher than the radiative ef- 
ficiency of such optically thin flows. For BZ mechanism, shocked flow produces higher 
efficiency than the shock free solutions and retrograde flow provides a slightly larger 
value of the efficiency than that for the prograde flow. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - shock 
waves - galaxies: active 
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1 INTRODUCTION 

Emission from the vicinity of black holes is usually con- 
sidered as powered by accretion. However, as proposed by 
Blandford & Znajek (1977), the kinetic energy of the spin- 
ning black hole may also provide another reservoir of the 
energy, which can be extracted with the help of the mag- 
netic field anchored in the surrounding plasma. Such energy 
is believed to power jets, and is expected to be small or 
comparable to the nominal accretion luminosity (Ghosh & 
Abramowicz 1997). 

The efficiency of the aforementioned Blandford-Znajek 
(BZ hereafter) mechanism is determined by the properties 
of the medium surrounding a black hole since the magnetic 
strength of the magnetic field and its respective topology re- 
flect the characteristic features of the surrounding medium. 

Most of the recent studies on BZ process have been per- 
formed for flow with considerably high angular momentum 
leading to the existence of a roughly Keplerian disk outside 
the Innermost Stable Circular Orbit (ISCO) and a plung- 
ing region closer to the horizon. Within this framework, BZ 
mechanisms is expected to be very efficient, although the de- 
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tails of the efficiency calculations differ in various works. Li 
& Paczynski (2000) proposed an alternating mechanism for 
disk accretion and black hole energy extraction which allows 
to enhance the pure accretion efficiency from 0.31 to 0.43. 
McKinney & Gammie (2004) in their numerical simulations 
of a disk with net dipole moment find efficiencies from 0.18 
up to 0.8. Garofalo (2009) suggests efficiency higher than 
0.0001 if the black hole spin is higher than 0.1, and reaching 
1 for classical solution, but for their flux-trapping scenario 
the corresponding numbers are shown to be much higher, 
0.01 and 10, respectively. 

Nemmen et al. (2007) and Garofalo (2009) addressed 
the issue of applicability of the BZ mechanism to nearby 
elliptical galaxies. However, the net angular momentum of 
such a quasi-Bondi inflow (as is believed to power such galax- 
ies) does not have to be high, and this can considerably 
change the maximum power available from the BZ energy 
extraction mechanism. 

Since the axisymmetric low angular momentum accre- 
tion onto a spinning black hole shows remarkably differ- 
ent features in comparison to a standard Keplerian/super- 
Keplerian disc accretion (see, e.g., Das & Czerny 2011, and 
references therein for further detail), there is a pressing need 
to study the energy extraction by BZ mechanism from low 
angular momentum rotating flow, both prograde and retro- 



2 Das & Czerny 



grade, to study the proper extraction mechanism of the spin 
energy of the massive black holes harboured by the nearby 
elliptical galaxies, and for our own Galactic centre black 
hole as well. Such mechanism may provide entirely differ- 
ent efficiency profile for the extracted energy, as well as for 
the corresponding luminosity function as compared with the 
other kind of galaxies where the accretion is of Keplerian in 
nature. Such difference may further finds its use as an im- 
portant tool to discriminate certain observational properties 
of various types of galaxies. 

Motivated by the aforementioned arguments, we present 
a detail analytical scheme of studying the efficiency of the 
BZ extracted energy from the general relativistic low angular 
momentum prograde and retrograde flow in the Kerr metric. 



2 LOW ANGULAR MOMENTUM 

ACCRETION WITH STATIONARY SHOCK 

The detail formulation of the equations governing the flow to 
be considered here and the procedure to obtain its general 
solutions are provided in Barai, Das & Wiita (2004) and 
Das & Czerny (2011), while the comprehensive characteristic 
features of such flow from the dynamical systems point of 
view are available in Goswami, Khan, Ray & Das (2006). 

The flow is assumed to be of negligible viscosity and 
the large radial inflow velocity, u, implies an almost free fall 
inflow close to the black hole and the advection effects are 
strong. Adiabatic accretion is considered. Two first integrals 
of motion, the conserved specific energy £ (the relativistic 
Bernoulli's constant) and the mass accretion rate M are ob- 
tained for such flow. The geometry of the flow is determined 
by a local flow thickness: function of the radial distance, 
sound speed and various other variables - that can be ob- 
tained by solving the general relativistic Euler equation in 
the vertical direction. The flow profile is essentially solved in 
the equatorial plane with vertically averaged values for the 
thermodynamic variables, and the complete stationary ac- 
cretion solution is specified by four astrophysicaly relevant 
parameters [£, A, 7, a], A, 7 and a being the specific flow an- 
gular momentum, the adiabatic index and the black hole 
spin, i.e., the Kerr parameter, respectively. 

The space gradient for the flow velocity u as well as 
the sound speed c s can be cast into a set of first order au- 
tonomous dynamical systems as: 
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from the differential solution of the energy conservation as 
well as the mass conservation equation of the following form: 
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where r, andp are the radial distance in the equatorial plane 
and the flow density. A and A are functions of r and a, vt is 
function of u, r, A, a, A, A, and other terms are functions of 
above quantities and of various metric elements as well as 
their space derivatives in the Boyer-Lindquist co-ordinate. 
Further details are available in Barai, Das & Wiita (2004) 
and Das & Czerny (2011). 

The fixed point analysis reveals that the flow described 
by the above model can have multiple (at most three) criti- 
cal points, with a centre type point flanked in between two 
saddle type points, for certain regions of the four dimen- 
sional parameter space spanned by [£,A,7,a]. Stationary 
shock may form if for such multi-critical flow the general 
relativistic Rankine Hugoniot condition gets satisfied. For 
flow with three critical points, the post shock flow will fi- 
nally pass through the inner sonic point while the the shock 
free solution passes only through the outer sonic point, as 
demonstrated in figure 2. of Das & Czerny (2011). The shock 
location, strength, and the ratio of various pre- to post shock 
dynamical and thermodynamic variables can be obtained as 
a function of the Kerr parameter, which indicates how the 
spin of the black hole influences such quantities. Dependence 
of such quantities on £,A and 7 can also be studied as well. 



3 CALCULATION OF EFFICIENCY FOR THE 
BZ MECHANISM 

We start with the general formula for the BZ luminosity 
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where £If is the angular velocity of the zero field observers, 
Qh is the angular velocity of the black hole horizon, rn is the 
horizon radius, Bh the normal component of the magnetic 
field. We calculate the efficiency of the energy extraction as 
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For flow variables calculated along the shocked flow 
passing through the inner sonic point, or along the shock 
free flow passing through the outer sonic point of the multi- 
transonic flow, or for flow through the only available ad- 
dle type sonic point for the mono-transonic flow, we de- 
fine the 'quasi-terminal value' (QTV hereafter) as Vs to be 
the value of any accretion variable V computed at a very 
close proximity rg to the event horizon, where r$ — r + + 5, 

r+ = 1 + (l- a 2 ) 3 and 8 = 0.001 GM J H , Mbh being the 
black hole mass. We integrate the flow upto r$ to obtain 
various QTV - Us, Ps,T$, Pg for example, where T and P 
are the flow temperature and pressure (both the gas pres- 
sure Pg as and total pressure P^°* = P| as + pgu 2 s can be 
"oomputed) respectively, and use such values to compute the 
i i^imum value of the efficiency factor r\ for the energy ex- 
traction corresponding to the BZ mechanism. 

We assume that the normal component of the magnetic 
field is in equipartition with the energy density since this 
maximizes the magnetic field supported by the inflowing 
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material. We also assume that £If = 0.5fin since it further 
maximizes the output (Phinney 1983, Reynolds et al. 2006). 
This assumption is not unrealistic, as shown by Komissarov 
(2001). 

In this work, quite often we compare the value of r\ for 
the shocked and shock free flow to understand the influence 
of shock on determining the value of 77. To accomplish that 
task, for same set of values of [£, A, 7, o], we calculate 77 from 
QTVs computed along the supersonic flow through the inner 
sonic point (77 for shocked flow) and for QTVs computed 
along the unperturbed supersonic flow through the outer 
sonic point (77 for shock free flow) and study the dependence 
of two such 77s on [£, A, 7, a] for prograde as well as for retro- 
grade flow. 

It is to be noted that in our work, for the first time in 
literature, low angular momentum accretion with shock has 
been studied for the counter rotating flow. 

Also to be noted that since for all values of [£, X, 7, a] 
the shock does not form, one requires to consider a specific 
limited range of fixed parameters to study the variation of 77 
on [£, A, 7, a] for flow with shock. All such issues are further 
clarified in section [4] using the respective figures. 

The value of 77 does not depend on the adopted value 
of the black hole mass and the QTV of the local density 
beyond the Bondi radius. However, it depends on the re- 
maining parameters: the black hole spin, asymptotic value 
of the Bernoulli constant, £, (or, equivalently, the tempera- 
ture), and the polytropic index, 7. 



4 RESULTS 

We calculate the stationary solutions for a low angular mo- 
mentum flow likely applicable to Sgr A* and black holes at 
the centers of elliptical galaxies. Specifically, we scale the 
solutions for the black hole mass and the asymptotic local 
density corresponding to the Sgr A*, 3.6 x 1O 6 M0, and 26 
cm~ 3 , respectively. For each of the solutions, the maximum 
efficiency of energy extraction through the BZ mechanism 
has been determined. 

In figure 1 we show the dependence of 77 on black hole 
spin a for prograde flow, for both shocked (red continuous 
line) and shock free (blue dashed line) solutions for fixed 
value of [£,7]. To cover a considerably large range for a 
(from slowly rotating to very fast rotation for the hole), we 
had to take three different values of A as mentioned in the 
figure, since flow characterized by a single value of A does 
not form the shock for the entire range of a analyzed here. 

It is obvious from the figure that the shocked flow pro- 
duces higher efficiency for energy extraction through the BZ 
mechanism since the post shock flow is relatively hotter and 
denser. This is a universal feature for dependence of 77 on any 
of initial boundary conditions corresponding to [£, A, 7, a], as 
is understood from the other figures as well. 

Similar results for the retrograde flow has been shown 
in figure 2. 

The dependence of 77 on various other flow parameters 
are shown in figure 3 for a fixed value of black hole spin as 
mentioned in the figure. One can have similar set of figures 
for any other value of a as well. The value of 77 correlates with 
£ and 7 since both factors lead to higher value of the pressure 
for a fixed density at the outer edge, which can be easily 
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Figure 1. The efficiency of energy extraction from the low an- 
gular momentum prograde accretion as a function of the black 
hole spin for solutions with (red continuous line) and without 
(blue dashed line) a shock for the angular momentum A equal 2.6 
(upper panel), 2.17 (middle panel) and 2.01 (lower panel). Fixed 
parameters are 7 = 1.43, £ = 1.00001. 



Figure 2. The efficiency of energy extraction from the low angu- 
lar momentum retrograde accretion as a function of the black hole 
spin, for solutions with (red continuous line) and without (blue 
dashed line) a shock for the A = 3.3, 7 = 1.4, and £ = 1.00001. 



seen in the case of Bondi flow. However, it anti-correlates 
with A for shocked case which is a manifestation of the fact 
that lower value of the angular momentum produces the 
shock closer to the event horizon and hence larger amount of 
energy can be extracted for such cases. 77, however, correlates 
with A for shock free transonic solutions. 

It is instructive to study whether the efficiency of the 
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energy extracted through the BZ mechanism is different 
for the counter rotating flow. One way to accomplish that 
task is to study the dependence of r\ on the same value 
(magnitude wise) of the Kerr parameter but for both the 
prograde and the retrograde flow, by keeping the values 
of [£,A,7] constant for those two cases. The findings has 
been presented in the figure 4 for a fixed set of values of 
[£ = 1.2, A = 2.0, 7 = 1.6]. Such high value of £ may not al- 
ways be realized for many astrophysical sources, however, 
such a value provides the transonic flow solution for the en- 
tire range of the Kerr parameter from a = — 1 to a = +1, 
and hence is useful for the aforesaid comparative study. Here 
the flow is mono-transonic and passes through a saddle type 
inner sonic point. 

The solutions show the monotonic decrease of the ra- 
dial position of the sonic point across the whole parameter 
range, — 1 < a < 1, changing from the value of 4.4 GMbh jc 2 
(a = —1.0) down to the horizon for a approaching 1. How- 
ever, other flow parameters, including pressure, achieve their 
minimum values at a — and increase with a for both pro- 
grade and retrograde flows. The rise, however, is not sym- 
metric due to the coupling with the sonic point asymmetry, 
as we illustrate in Fig. [4] The flow efficiency is almost sym- 
metric, with the values 7 x 10 -3 for a = -0.99 and 6 x 10~ 3 
for a — 0.99, since the dominant factor is the a 2 dependence, 
as the flow density close to the horizon is proportional to the 
pressure, which simply reflects the fact that the sound speed 
reaches quite close to c at the horizon. 



5 DISCUSSION 

The extraction of the rotational energy from the black hole 
is an interesting option to power jet outflows from direct 
vicinity of a black hole (Blandford & Znajek 1977). The 
efficiency of this mechanism is currently under thorough in- 
vestigation, and it depends on the Kerr parameter as well as 
the properties of the inflowing material and the surrounding 
large scale magnetic field. 

Most of the attention was payed to the matter in disk- 
like accretion since it this case the efficiency may be the 
maximum. As shown by Tchekhovskoy et al. (2011), this 
efficiency may considerably be higher than the accretion ef- 
ficiency, thus exceeding 100 per cent. 

However, there are observational arguments that jets 
preferentially form when the accretion pattern is different 
from the classical Shakura-Sunyaev optically thick, geomet- 
rically thin flow (e.g. Fender et al. 2004). Such flows are 
generally known as RIAF (Radiatively Inefficient Accretion 
Flows) and they can be divided into high accretion angu- 
lar momentum flows (usually known as ADAF models) and 
low angular momentum accretion flows, including spheri- 
cally symmetric Bondi accretion. 

In this paper we analysed for the first time the effi- 
ciency of the Blanford-Znajek mechanism in the case of low 
angular momentum flows with, and without shocks which 
may be applicable to elliptical galaxies or to Sgr A*. For 
the first ever time in the literature, we have also studied the 
shocked flow and the related efficiency of BZ mechanism for 
the retrograde flow as well. This is important because the 
relevance of the counter rotating flow in black hole astro- 
physics is being increasingly evident from recent theoretical 



O 




5xl0" 6 
log 8 



tafj 
O 



QTj 
O 



-2.8 
-3 
-3.2 
-3.4 

-2.8 
-3 
-3.2 
-3.4 



1 — i — | — i — i — i — i — | — i — i — i — i — | — i — r 



j i I i i i i I i i i i I i L 



1.41 1.42 1.43 

7 



I — i — i — r 



I — i — i — r 



j i i i_ 



j i i L 



_L 



j i i i_ 



2.15 



2.2 



2.25 



2.3 



A 



Figure 3. The efficiency of energy extraction from the low an- 
gular momentum accretion flow as a function of the Bernoulli 
constant (upper panel), polytropic index (middle panel) and an- 
gular momentum (lower panel) for solutions with (red continuous 
line) and without (blue dashed line) shock. Fixed parameters are 
7 = 1.43, A = 2.17 (upper panel), £ = 1.000004, A = 2.29 (middle 
panel),£ = 1.000004, 7 = 1.43 (lower panel), a = 0.881049812 (all 
panels). 



and observational findings (Nixon et al. 2010, Dauser et al. 
2011 and references therein). 

We find the efficiency of the energy extraction in this 
case to be of the order of 0.1 %, generally somewhat lower 
that in ADAF type of inflow, and much lower than in the 
thin disk accretion. Narayan & Fabian (2011) considered a 
model with an inflow timescale only a few times the free- 
fall time, due to a combination of a moderate initial angular 
momentum of the material and very high viscosity, and they 
obtained energy extraction efficiency of order of up to 2%. 
However, the expected radiative efficiency of the flow con- 
sidered in our paper is also very low, at the level of 10 -6 , so 
the BZ energy extraction can easily dominate in such flows, 
as for low accretion rate ADAF solutions (e.g. Wu & Cao 
2008). 

There is also an additional mechanism which can in- 
crease the energy extraction efficiency. Narayan et al. (2003) 
postulated that the magnetic field can accumulate close to 
the black hole horizon. In this case the support for the mag- 
netic filed can come also from the ram pressure, and from 
the maximum pressure at the flow, as discussed by Narayan 
& Fabian (2011). If we include the ram pressure into our flow 
model, the efficiency is enhanced by a factor ~ 15 for pro- 



Blandford- Znajek mechanism for low angular momentum accretion 5 



i — i — i — i — I — i — i — i — r 



l — i — i — I — i — i — i — r 



6.5 - 



(X 

g> 6 



+ 

T 
I- 

I 

I _ 
/ - 



5.5 




grade shocked flow, as shown in Fig. [5] (similar results can be 
obtained for the retrograde flow as well), but it never reaches 
the efficiency higher than 100 % seen by Tchekhovskoy et 
al. (2011) due to relatively low angular momentum assumed 
in our model. 

Our consideration of the energy extraction in low an- 
gular momentum flow shows a positive correlation with the 
spin, and the efficiency for a = 0.99 is by a factor 10 higher 
than for a = 0.5. Similar trend is seen in the flows based on 
other description of the dynamics. Some observations, how- 
ever, do not confirm any obvious correlation between the jet 
power and the measured black hole spin in active galaxies 
(Broderick & Fender 2011) or black hole binaries (Fender et 
al. 2010). 

The role of the BZ energy extraction mechanism in jet 
formation is thus still an open question, and a combination 
of two mechanisms is quite likely, as concluded by Foschini 
(2011), with an increasing role of rotational energy extrac- 
tion with decreasing object luminosity (Armitage & Natara- 
jan 1999), including Bondi and low angular momentum flow. 
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Figure 4. The pressure close to the horizon and the efficiency 
of energy extraction from the low angular momentum accretion 
flow for solutions without a shock covering the whole range of co- 
rotating and counter-rotating black holes. Fixed parameters are 
A = 2.0, 7 = 1.6, 8 = 1.2. 




a 

Figure 5. The efficiency of the energy extraction when the ram 
pressure is included (blue long dashed line) in comparison with 
the gas pressure only (red continuous line), [£, A, 7] used are same 
as in the panel one of Fig. 1. 
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ABSTRACT 

Blandford-Znajek (BZ) mechanism has usually been studied in the literature for ac- 
cretion with considerably high angular momentum leading cither to the formation 
of a cold Keplerian disc, or a hot and geometrically thick sub-Keplerian flow as de- 
scribed within the framework of ADAF/RIAF. However, in nearby elliptical galaxies, 
as well as for our own Galactic centre, accretion with very low angular momentum 
is prevalent. Such quasi-spherical strongly sub-Keplerian accretion has complex dy- 
namical features and can accommodate stationary shocks. In this letter, we present 
our calculation for the maximum efficiency obtainable through the BZ mechanism for 
complete general relativistic weakly rotating axisymmetric flow in the Kerr metric. 
Both shocked and shock free flow has been studied in detail for rotating and counter 
rotating accretion. Such study has never been done in the literature before. We find 
that the energy extraction efficiency is low, about 0.1%, and increases by a factor 15 
if the ram pressure is included. Such an efficiency is still much higher than the radia- 
tive efficiency of such optically thin flows. For BZ mechanism, shocked flow produces 
higher efficiency than the shock free solutions and retrograde flow provides a slightly 
larger value of the efficiency than that for the prograde flow. 
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1 INTRODUCTION 

Emission from the vicinity of black holes is usually con- 
sidered as powered by accretion. However, as proposed by 
Blandford & Znajek (1977), the kinetic energy of the spin- 
ning black hole may also provide another reservoir of the 
energy, which can be extracted with the help of the mag- 
netic field anchored in the surrounding plasma. Such energy 
is believed to power jets, and is expected to be small or 
comparable to the nominal accretion luminosity (Ghosh & 
Abramowicz 1997). 

The efficiency of the aforementioned Blandford-Znajek 
(BZ hereafter) mechanism is determined by the properties 
of the medium surrounding a black hole since the magnetic 
strength of the magnetic field and its respective topology re- 
flects the characteristic features of the surrounding medium. 

Most of the recent studies on BZ process have been per- 
formed for flow with considerably high angular momentum 
corresponding to the circularization radius of order of hun- 
dreds gravitational radii. Accretion is then possible due to 
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viscosity, and the angular momentum redistribution leads 
to the existence of a cold Keplerian or a hot sub-Keplerian 
disk outside the Innermost Stable Circular Orbit (ISCO), 
and a plunging region closer to the horizon. Within this 
framework, BZ mechanisms is expected to be very efficient, 
although the details of the efficiency calculations differ in 
various works. Li & Paczynski (2000) proposed an alternat- 
ing mechanism for cold disk accretion and black hole energy 
extraction which allows to enhance the pure accretion ef- 
ficiency from 0.31 to 0.43. Hot flows are also likely to be 
very efficient in energy extraction (e.g. Meier 2001, De Vil- 
liers et al. 2003, McKinney 2006, Hawley & Krolik 2006, 
Nemmen et al. 2007, Wu & Cao 2008, Tchekhovskoy et al. 
2010). McKinney & Gammie (2004) in their numerical sim- 
ulations of a disk with net dipole moment find efficiencies 
up to 0.068. Garofalo (2009) suggests efficiency higher than 
0.0001 if the black hole spin is higher than 0.1, and reaching 
1 for classical solution, but for their flux-trapping scenario 
the corresponding numbers are shown to be much higher, 
0.01 and 10, respectively. 

Nemmen et al. (2007) and Garofalo (2009) addressed 
the issue of applicability of the BZ mechanism to nearby 
elliptical galaxies. However, the net angular momentum of 
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such a quasi-Bondi inflow (as is believed to power such galax- 
ies) does not have to be high. The gravitationally relaxed 
isotropic distribution of the old population stars at the cen- 
tre cluster of the Milky way (Ott et al. 2003), for example, 
may contribute the bulk of the infalling material, leading to 
the low angular momentum accretion flow. Currently, there 
is no reliable estimates of the angular momentum of the ma- 
terial reaching the vicinity of the black hole in such systems. 
The best simulations performed for Sgr A* (Cuadra et al. 
2008) imply quite a large circularization radius of 5000 grav- 
itational radii but calculated at the inner edge of the com- 
putational domain, and most of the modeled clumpy gas is 
unlikely to reach the black hole. Thus the net angular mo- 
mentum may be much lower. This scenario can considerably 
change the maximum power available from the BZ energy 
extraction mechanism. 

Since the axisymmetric low angular momentum accre- 
tion onto a spinning black hole shows remarkably differ- 
ent features in comparison to a standard Keplerian/super- 
Keplerian disc accretion (see, e.g., Das & Czerny 2012, and 
references therein for further detail) , there is a pressing need 
to study the energy extraction by BZ mechanism from low 
angular momentum rotating flow, both prograde and retro 
grade, to study the proper extraction mechanism of the spin 
energy of the massive black holes harboured by the nearby 
elliptical galaxies, and for our own Galactic centre black 
hole as well. Such mechanism may provide entirely differ- 
ent efficiency profile for the extracted energy, as well as of 
the corresponding luminosity function as compared with the 
other kind of galaxies where the accretion is of Keplerian in 
nature. Such difference may further finds its use as an im- 
portant tool to discriminate certain observational properties 
of various types of galaxies. 

Motivated by the aforementioned arguments, we present 
a detail analytical scheme of studying the efficiency of the 
BZ extracted energy from the general relativistic low angular 
momentum prograde and retrograde flow in the Kerr metric. 



2 LOW ANGULAR MOMENTUM 

ACCRETION WITH STATIONARY SHOCK 

Our general relativistic low angular momentum accre- 
tion possesses dynamical features intermediate between the 
Bondi (Bondi 1952) type accretion and ADAF (Advection 
Dominated Accretion Flow, see, e.g., Narayan & Yi 1994). 
The angular momentum at the outer flow boundary is zero 
in spherical flow, is locally a fraction of the Keplerian an- 
gular momentum in ADAF, and it is of order of the Kep- 
lerian angular momentum at ISCO for the flow model con- 
sidered here. The corresponding circularization radius is not 
too large in comparison to the innermost stable circular or- 
bit (ISCO), hence the salient features of such flow in quite 
similar to the spherical flow at the outer region whereas sig- 
nificant departure in the inner region is realized through the 
appearance of the angular momentum barrier and station- 
ary shock may form as a consequence. Thus our model, as 
compared to Bondi flow as well as ADAF/RIAF, however, 
has richer dynamical features, allowing the flow to become 
multi-transonic and to accommodate a standing shock. 

The detailed formulation of the equations governing the 
flow to be considered here and the procedure to obtain its 



general solutions analytically are provided in Barai, Das & 
Wiita (2004) and Das & Czerny (2012), while the compre- 
hensive characteristic features of such flow from the dynam- 
ical systems point of view are available in Goswami, Khan, 
Ray & Das (2006). The numerical HD and MHD simula- 
tions of such a flow were performed by Proga & Begel- 
man (2003a, b). Here we use analytical solutions to cover 
the whole parameter range. 

The flow is assumed to posses considerable advective 
(radial three- velocity) velocity u to make the infall time scale 
smaller than the viscous time scale in order to validate the 
inviscid approximation. Adiabatic accretion is considered. 
Two first integrals of motion, the conserved specific energy 
£ (the relativistic Bernoulli's constant) and the mass accre- 
tion rate M are obtained for such flow. The geometry of the 
flow is determined by a local flow thickness: function of the 
radial distance, sound speed and various other variables - 
that can be obtained by solving the general relativistic Euler 
equation in the vertical direction. The flow profile is essen- 
tially solved on the equatorial plane with vertically averaged 
values for the thermodynamic variables, and the complete 
stationary accretion solution is specified by four astrophys- 
icaly relevant parameters [£,X, 7, a], A, 7 and a being the 
specific flow angular momentum, the adiabatic index and 
the black hole spin, i.e., the Kerr parameter, respectively. 

The space gradient for the flow velocity u as well as 
the sound speed c s can be cast into a set of first order au- 
tonomous dynamical systems as: 
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from the differential solution of the energy conservation as 
well as the mass conservation equation of the following form: 
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H being the local flow thickness: 
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where r, andp are the radial distance on the equatorial plane 
and the flow density. A and A are functions of r and a, vt is 
function of u, r, A, a, A, A, and other terms are functions of 
above quantities and of various metric elements as well as 
their space derivatives in the Boyer-Lindquist co-ordinate. 
Further details are available in Barai, Das & Wiita (2004) 
and Das & Czerny (2012). 

The fixed point analysis reveals that the flow described 
by the above model can have multiple (at most three) criti- 
cal points, with a centre type point flanked in between two 
saddle type points, for certain regions of the four dimen- 
sional parameter space spanned by [£, A, 7, a]. Stationary 
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shock may form if for such multi-critical flow the general 
relativistic Rankine Hugoniot condition gets satisfied. For 
flow with three critical points, the post shock flow will fi- 
nally pass through the inner sonic point while the the shock 
free solution passes only through the outer sonic point, as 
demonstrated in figure 2. of Das & Czerny (2012). The shock 
location, strength, and the ratio of various pre- to post shock 
dynamical and thermodynamic variables can be obtained as 
a function of the Kerr parameter, which indicates how the 
spin of the black hole influences such quantities. Dependence 
of such quantities on £,A and 7 can also be studied as well. 



-3 



3 CALCULATION OF EFFICIENCY FOR THE 
BZ MECHANISM 

We start with the general formula for the BZ luminosity 

1 Qf(SIh — Of) d 2 2 2 ,„s 
= 32 fi2 B nr H a c, (6) 

where £1f is the angular velocity of the zero field observers, 
Qh is the angular velocity of the black hole horizon, rn is the 
horizon radius, Bh the normal component of the magnetic 
field. We calculate the efficiency of the energy extraction as 



For flow variables calculated along the shocked flow 
passing through the inner sonic point, or along the shock 
free flow passing through the outer sonic point of the multi- 
transonic flow, or for flow through the only available ad- 
dle type sonic point for the mono-transonic flow, we de- 
fine the 'quasi-terminal value' (QTV hereafter) as Vs to be 
the value of any accretion variable V computed at a very 
close proximity rg to the event horizon where r& — r+ + 5, 

r+ = 1 + (1 - a 2 ) ^ and S = 0.001 GM J H , Mbh being the 
black hole mass to be considered. We integrate the flow upto 
r$ to obtain various QTV - Us, ps,Ts, Ps for example, where 
T and P are the flow temperature and pressure (both the 
gas pressure Pg a3 and total pressure Pg ot = Pg a3 + psiig can 
be computed) respectively, and use such values to compute 
the maximum value of the efficiency factor 77 for the energy 
extraction corresponding to the BZ mechanism. 

We assume that the normal component of the magnetic 
field is in equipartition with the energy density since this 
maximizes the magnetic field supported by the inflowing 
material. We also assume that VLf — 0.5Qh since it further 
maximizes the output (Phinney 1983, Reynolds et al. 2006). 
This assumption is not unrealistic, as shown by Komissarov 
(2001). 

In this work, we frequently compare the value of 77 for 
the shocked and shock free flow to understand the influence 
of shock on determining the value of 77. To accomplish that 
task, for same set of values of [£, A, 7, o], we calculate 77 from 
QTVs computed along the supersonic flow through the inner 
sonic point (77 for shocked flow) and for QTVs computed 
along the unperturbed supersonic flow through the outer 
sonic point (77 for shock free flow) and study the dependence 
of two such 77s on [£, A, 7, a] for prograde as well as for retro- 
grade flow. 
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Figure 1. The efficiency of energy extraction from the low angu- 
lar momentum prograde accretion as a function of the black hole 
spin for solutions with (red continuous line) and without (blue 
dashed line) a shock. For the entire region of the black hole spin 
a shown in the figure, the angular momentum A varies from 2.6 
(lower values of a) down to 2.01 (higher values of a) to ensure the 
existence of stationary solutions. Fixed parameters are 7 = 1.43, 
£ = 1.00001. 
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Figure 2. The efficiency of energy extraction from the low angu- 
lar momentum retrograde accretion as a function of the black hole 
spin, for solutions with (red continuous line) and without (blue 
dashed line) a shock for the A = 3.3, 7 = 1.4, and £ = 1.00001. 



It is to be noted that in our work, for the first time in 
literature, low angular momentum accretion with shock has 
been studied for the counter rotating flow. 

Since the stationary solutions exist only for a subset of 
values of [£, A, 7, a], and the solutions with shocks cover even 
smaller parameter space, one requires to consider a specific 
limited range of fixed parameters to study the variation of 77 
on [£, A, 7, a] for flow with shock. All such issues are further 
clarified in section [4] using the respective figures. 

The value of 77 does not depend on the adopted value 
of the black hole mass and the QTV of the local density 
beyond the Bondi radius. However, it depends on the re- 
maining parameters: the black hole spin, asymptotic value 
of the Bernoulli constant, £, (or, equivalently, the tempera- 
ture), and the polytropic index, 7. 
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4 RESULTS 

We calculate the stationary solutions for a low angular mo- 
mentum flow likely applicable to Sgr A* and black holes at 
the centers of elliptical galaxies. Specifically, we scale the 
solutions for the black hole mass and the asymptotic local 
density corresponding to the Sgr A*, 3.6 x 1O 6 M0 (Eckart et 
al. 2006), and 26 cm -3 (Baganoff et al. 2003), respectively. 
For each of the solutions, the maximum efficiency of energy 
extraction through the BZ mechanism has been determined. 

In figure I we show the dependence of rj on black hole 
spin a for prograde flow, for both shocked (red continuous 
line) and shock free (blue dashed line) solutions for fixed 
value of [£,"/]■ To cover a considerably large range for a 
(from slowly rotating to very fast rotation for the hole), we 
had to take a range of values of A as mentioned in the figure 
caption, since flow characterized by a single value of A does 
not form the shock for the entire range of a analyzed here. 

It is obvious from the figure that the shocked flow pro- 
duces higher efficiency for energy extraction through the BZ 
mechanism since the post shock flow is relatively hotter and 
denser. This is a universal feature for dependence of rj on any 
of initial boundary conditions corresponding to [£, A, 7, a], as 
is understood from the other figures as well. Similar results 
for the retro grade flow has been shown in figure 2. 

The dependence of rj on £, A and 7 is rather weak - it 
would create a finite but small width to the plot shown in 
Figs. 1 and 2. For a given value of a, varying A gives the 
change in log 77 by less than 0.1, varying 7 gives a change 
less than 0.2, and varying £ from 1. 000001 to 1.00001 gives 
a change less than 0.1. Much larger values of £ can increase 
the efficiency but it imply unrealistically high temperature 
at large distances from the black hole. The value of 77 cor- 
relates with £ and 7 since both factors lead to higher value 
of the pressure for a fixed density at the outer edge, which 
can be easily seen in the case of Bondi flow. However, it 
anti-correlates with A for shocked case which is a manifesta- 
tion of the fact that lower value of the angular momentum 
produces the shock closer to the event horizon and hence 
larger amount of energy can be extracted for such cases, rj, 
however, correlates with A for shock free transonic solutions. 

It is instructive to study whether the efficiency of the 
energy extracted through the BZ mechanism is different for 
the counter rotating flow. One way to accomplish that task 
is to study the dependence of rj on the same value (mag- 
nitude wise) of the Kerr parameter but for both the pro- 
grade and the retrograde flow, by keeping the values of 
[£,A,7] constant for those two cases. Such observation has 
been presented in the figure 3 for a fixed set of values of 
[£ — 1.2, A = 2.0, 7 = 1.6]. Such high value of £ may not al- 
ways be realized for many astrophysical sources, however, 
such a value provides the transonic flow solution for the en- 
tire range of the Kerr parameter from a — —1 to a = +1, 
and hence is useful for the aforesaid comparative study. Here 
the flow is mono-transonic and passes through a saddle type 
inner sonic point. 

The solutions show the monotonic decrease of the ra- 
dial position of the sonic point across the whole parameter 
range, —1 < a < 1, changing from the value of 4.4 GMbh/c 2 
(a = —1.0) down to the horizon for a approaching 1. How- 
ever, other flow parameters, including pressure, achieve their 
minimum values at a — and increasing with a for both pro- 




Figure 3. The pressure to the density ratio close to the hori- 
zon and the efficiency of energy extraction from the low angular 
momentum accretion flow for solutions without a shock covering 
the whole range of co-rotating and counter-rotating black holes. 
Fixed parameters are A = 2.0, 7 = 1.6, £ = 1.2. Dotted line 
marks the a 2 dependence. 



grade and retrograde flows. The rise, however, is not sym- 
metric due to the coupling with the sonic point asymmetry, 
as we illustrate in Fig. [3] The flow efficiency is almost sym- 
metric, with the values 7 x 10~ 3 for a = -0.99 and 6 x 10~ 3 
for a = 0.99, since the dominant factor is the a 2 dependence, 
as shown using the dotted line, as the flow density close to 
the horizon is proportional to the pressure, which simply re- 
flects the fact that the sound speed reaches quite close to c 
at the horizon. 



5 DISCUSSION 

The extraction of the rotational energy from the black hole 
is an interesting option to power jet outflows from direct 
vicinity of a black hole (Blandford & Znajek 1977). The 
efficiency of this mechanism is currently under thorough in- 
vestigation, and it depends on the Kerr parameter as well as 
the properties of the inflowing material and the surrounding 
large scale magnetic field. 

Most of the attention was payed to the matter in disk- 
like accretion since it this case the efficiency may be the 
maximum. As shown by Tchekhovskoy et al. (2011), in the 
extreme case of a magnetically-arrested flow this efficiency 
may considerably be higher than the accretion efficiency, 
thus exceeding 100 per cent. 

However, there are observational arguments that jets 
preferentially form when the accretion pattern is different 
from the classical Shakura-Sunyaev optically thick, geomet- 
rically thin flow (e.g. Fender et al. 2004). Such flows are 
generally known as RIAF (Radiatively Inefficient Accretion 
Flows; see Quataert 2001) and they can be divided into high 
accretion angular momentum flows (usually known as ADAF 
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models) and low angular momentum accretion flows, includ- 
ing spherically symmetric Bondi accretion. 

In this paper we analyzed for the first time the efficiency 
of the Blandford-Znajek mechanism in the case of low angu- 
lar momentum flows. Such inflow starts at the outer bound- 
ary with an angular momentum much smaller than the lo- 
cal Keplerian value and proceed almost as a free fall till the 
angular momentum barrier becomes important. Stationary 
solutions with negligible viscosity exist for a range of param- 
eters, with and without shocks. These solutions may be ap- 
plicable to elliptical galaxies or to Sgr A*. For the first ever 
time in the literature, we have also studied the BZ mech- 
anism for a shocked flow, and the low angular momentum 
solutions for the retrograde flow. This is important because 
the relevance of the counter rotating flow in black hole astro- 
physics is being increasingly evident from recent theoretical 
and observational findings (Nixon et al. 2010, Dauser et al. 
2011 and references therein). The retrograde flow is more 
efficient by ~ 10 %, so the enhancement is much lower that 
in the model considered by Garofalo et al. (2007). 

We find the efficiency of the energy extraction to be of 
the order of 0.1 %, generally somewhat lower that in ADAF 
type of inflow, and much lower than in the thin disk accre- 
tion. Narayan & Fabian (2011) considered a model with an 
inflow timescale only a few times the free-fall time, due to a 
combination of a moderate initial angular momentum of the 
material and very high viscosity, and they obtained energy 
extraction efficiency of order of up to 2%. In general, for a 
given accretion rate, the inflow velocity in the outer part of 
ADAF is by a factor a lower than the free fall, where a is the 
viscosity parameter, and the density is by a factor 1/a higher 
than in low angular momentum flow in the outer parts, and 
this difference is only partially diminished in the innermost 
part, particularly in the presence of shocks. However, the 
expected radiative efficiency of the flow considered in our 
paper is also very low, emissivity calculated for a particular 
ex ample of a transonic flow studied in the post Newtonian 
framework was of order of 7 x 10~ 6 (Moscibrodzka et al. 
2006), and hence the BZ energy extraction can easily domi- 
nate in such flows, as for low accretion rate ADAF solutions 
(e.g. Wu & Cao 2008). The exact radiative efficiency both in 
ADAF and low angular momentum flow depends crucially 
on assumptions about the electron-ion coupling, but the ra- 
diative emissivity (synchrotron and bremsstrahlung) scale 
roughly quadratically with density under the same assump- 
tions so the radiative efficiency of low angular momentum 
is lower roughly by a factor 1/a 2 . Thus in the low angu- 
lar momentum flow the BZ energy extraction more easily 
dominates the radiative efficiency than in ADAF. 

Our low angular momentum model does not include the 
effect of the viscosity since the flow possesses considerable 
radial advection, for which the viscous time scale is fairly 
large compared to the infall time scale. The stationary so- 
lutions used in this paper would not be much affected if 
the viscosity is included. On the other hand, large viscos- 
ity would likely broaden the parameter range for stationary 
solutions and the model would become similar to ADAF. 

There is also an additional mechanism which can in- 
crease the energy extraction efficiency. Narayan et al. (2003) 
postulated that the magnetic field can accumulate close to 
the black hole horizon. In this case the support for the mag- 
netic filed can come also from the ram pressure, and from 



the maximum pressure at the flow, as discussed by Narayan 
& Fabian (2011). If we include the ram pressure into our 
flow model, the efficiency is enhanced by a factor ~ 15, and 
the model reaches the jet efficiency of ~ 2 % as required 
by the Chandra observations of the elliptical galaxies (Allen 
et al. 2006, Nemmen et al. 2007), if the black hole spin is 
larger than 0.9 for realistic value of the Bernoulli constant 
(Fig. 1). Thus with the ram pressure included, the BZ mech- 
anism works equally efficiently in low angular momentum 
model as in ADAF. However, it never reaches the efficiency 
higher than 100 % seen by Tchekhovskoy et al. (2011) due 
to accumulation of the magnetic field close to the horizon, 
so the magnetic field pressure dominates by two orders of 
magnitude the ram pressure in their model. 

Our consideration of the energy extraction in low angu- 
lar momentum flow shows a dependence on the spin some- 
what steeper than the simple a 2 law, but the enhancement 
in comparison with this law is only by a factor 2 - 3 at the 
a — 1 limit. It is less than the factor ~ 100 in the numerical 
results of McKinney (2005) for the BZ jet efficiency from a 
thick disk. It may partially reflect the difference in the dy- 
namics, but our results may slightly underestimate the effect 
due to assumption of the field equipartition while McKinney 
(2005) and McKinney & Gammie (2004) compute the mag- 
netic field self-consistently. Some observations, however, do 
not confirm any obvious co relation between the jet power 
and the measured black hole spin in active galaxies (Brod- 
erick & Fender 2011) or black hole binaries (Fender et al. 
2010). 

The role of the BZ energy extraction mechanism in jet 
formation is thus still an open question, and a combination 
of two mechanisms is quite likely, as concluded by Foschini 
(2011), with an increasing role of rotational energy extrac- 
tion with decreasing object luminosity (Armitage & Natara- 
jan 1999), including Bondi and low angular momentum flow. 
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